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^ ' The static current operator leads to definitional zero frequency divergence and unphysical results 

, in studying nonlinear optical susceptibilities of polymers. A well-defined dipole-dipole correlation 

' is superior to the complicated current-current correlation to solve this problem. As illustrative 

I examples, optical susceptibilities under both SSH and TLM models of trans-(CH)i, are studied. 

■ New analytical results are obtained. The reasons of previous improper results are analyzed. 
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-4-^ ^ To study the nonlinear optical (NLO) properties of polymers, periodic approximate models are necessary to simplify 
' the real systems. Some good approximate models of polymers based on the tight-binding approximation (TBA), such 
"j^ [ as SSH 1^ and TLM |^] models in weakly correlated systems, Hubbard and Pariser-Parr-Pople models in strongly 
^ ■ correlated and electron-electron (e-e) interaction systems, etc, have yielded physical insights surpassing the compli- 
I ' cated non-approximate computations. In considering the optical response of these models, a U(l) transformation has 
■ O been suggested to provide the gauge invariance of the TBA Hamiltonian. 

^ ' In linear optical (LO) response theory, the Kubo formula based on current-current (J-J) correlation is widely 
^ used because of simplicity. It is commonly held that in discussing the optical susceptibilities of materials, the J-J 
I— 1^ correlation (p • A) will play the same role as that of dipole-dipole (D-D) correlation (E • r) |^,^ and that the apparent 
^ ■ definitional zero frequency divergence (ZFD) in the NLO susceptibilities definition |10jl^ is only a virtual problem, 
^ ' although the proofs only have been shown under certain assumptions. However, there exists some discrepancies 

\^ from the above conclusion under the models: (i) In the study of conductivity a under TBA models, |13jlJ] the J-J 
correlation will have ZFD in Im(cr) (related to the real part of J-J) although it can give the correct Re(a) (related 
\l to the imaginary part of J-J). (ii) In third-harmonic generation (THG) of trans-(CH)a;, the experimentally observed 
two-photon absorption peak (TP) ||l^ has given rise to the theoretical explanations. Based on J-J correlation, the 
TP was explained [lllj by choosing the TLM model. But it has been criticized by others [l^-|2C|] with simple parity 
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■ consideration. Also the spectrum is quite different from that under the D-D correlation. Thus, it casts some 
"Y^ ■ doubts in the practical application of J-J correlation. 

Ci ' In this letter, we will show that all above controversies are caused by the improper application of J-J correlation 
^ ■ based only on the static current operator Jq. To recover the correct results, a well-defined D-D correlation is more 
' ' suitable for studying the optical susceptibilities than the J-J correlation, whose equivalence to the D-D correlation 
^ . can be satisfied by introducing complicated induced field currents (IFCs). Besides solving the apparent ZFD in the 
Q ' definition of NLO susceptibility, the D-D correlation is to be favored over the J-J correlation due to the lack of the 
O ■ gauge dependence of the vector potential A |^ and simplifying the definitional complexity of the contribution by 
IFCs in the high order expansions. j3|Q| As a deduction, the TP cusp and ZFD ||ll|,|l^,Q will no longer exist in both 
SSH and TLM models by the well-defined D-D correlation. However, this seemingly trivial conclusion has not been 
clearly illustrated by the others although the D-D correlation has already been applied in the NLO response of the 
^ : real systems, g^,^ 

" " ' The application of J-J correlation is simplier than that of the D-D correlation in LO response because of the 
convenient search for the static current Jq compared with that for dipole expression D under the approximate 
models. As a typical example, the position operator r is ill-defined in periodic systems p^,|2^ i n real space while Jq 
is not. Further, the static current Jq still can give the correct results in the LO absorption |7| J13| , p^ because the IFCs 
only contribute to the real part of J-J correlation. Thus the simplicity of the Kubo formula is still satisfied by applying 
the static current Jo and ZFD in the real part of J-J correlation could be solved either by the Kramers-Kronig (KK) 
relation or by subtracting {[j,j]){uj = 0). [|l^ It is not a big surprise to see why the correct results still can be 
preserved for the LO response through Jq. However, the advantage of J-J correlation application in LO response is 
no longer available for the NLO response because of the IFCs, 0,|l3| thus the definition of n-th order J-J correlation 
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based only on the static current ||Tl[ will lead to the incorrect results. Besides the difficulty in obtaining the correct 
IFCs in approximate models, J- J correlation will make the definition of NLO susceptibilities unable to be written in 
a general form as defined in Wu' work and will make the computations very tedious, although it may give the 
correct result. 

In periodic systems, the imposed boundary condition requires (i) Bloch wave functions to describe the extended 
electronic states and (ii) all operators including the position operator f to have at least the same periodicity. ||2|] A 
band index n and crystal momentum k are used to label the Bloch states: |n,k >= u„.k(r)e'''''', where Un,k(r) is the 
periodic function under the translation of lattice vector. If the homogeneous electric field is applied, the energy level 
En{k) should be changed to En{n) [pf and the static current operator Jo(k) should be replaced by a new current 
operator J(k), where k = k — eA/h. Iiq] By Taylor expansion of J(k) based on powers of the vector potential A, we 
will obtain the IFCs related to all orders of A besides the static current Jq. There is no doubt that the IFCs will 
contribute to the optical response. Complexities in handling IFCs (usually related to the intra-band currents |jl^) in 
NLO studies could be imagined. Polymer models usually share the same properties as periodic systems. The 
direct consequence of dropping the IFCs will result in ZFD and unphysical results. ||ll|,|l3|,0 Although the unphysical 
results in NLO response have been recurrently questioned by the others, the direct reason has never been 

revealed and correct analytical results have never been obtained. The chief difficulty is that direct ways to include IFCs 
in many approximate models are forbidden and the static current could easily be improperly used. By 

a well-defined D-D correlation, we will show that all difficulties caused by the improper use of J- J correlation could 
be solved. To give an intuitive picture, we will illustrate those correlations by studying the simple one-dimensional 
TBA electron-lattice models of trans-(CII)a; (SSH and TLM) whose optical properties have been widely studied 
by others. ||ll| , [l^ - pO|,p^ To avoid the ill-definition and to provide the periodicity of the position operator f , we express 



f under |n, k > as: 

r„k,«'k' = i5n,n'^kSik - k') + f7„,„,(k)5(k - k'), (1) 



i;.(r)VkW„'^k(r)(ir and v is unit cell volume. 



Susceptibilities definition by D-D correlation: Without considering retardation effect, ^ optical susceptibilities 
are usually defined by expanding the optical polarization in powers of the transverse electric field. Under that 
approximation, the general nth-order susceptibility is a purely material quantity defined as: M 
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Jdr,-- dr„ Jdt,--- dt^ J dvdt e"*"-+*^"(f D(r, t)t){r^M) ■ ■ ■ D(r„, i„)), (2) 



where = — w;, T is the time-ordering operator and D is dipole operator. 

■i=i 

Based on periodic TBA, SSH Hamiltonian M IS given: 



HssH — - ^ 



io + (-l)'f' 



where ta is the transfer integral between the nearest-neighbor sites, A is the gap parameter and Cj ^{Ci^s) cre- 
ates (annihilates) an vr electron at site I with spin s. In contimmm limitation, above SSH model will give the TLM 
model. II 

In momentum space, the above Hamiltonian with electron-photon (e-A) interaction could be found as follows: 

H{k, t)^Yl £ik)i'ijt)a3i^kAt) - D ■ Ee'^\ (3) 



where e{k) — y [2tQC0s{ka)]^ + [Asin{ka)]'^ and -0^ si^)~ i^k'^si^)^ '^c'sW) two-component spinor describing ex- 

citations of electrons in the conduction band and valence band. Long wave approximation 0| is applied in electro- 
magnetic field E with frequency uj. The dipole operator D could be obtained by the Eq. (|l|): 

i) = e^(/3(fc) i^l,a2i^k.s + ^^ i^ljk.s), (4) 

k.s 
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where /?(fc) = — Atoa/e^(fc), is the coefficient related to the interband transition between the conduction and valence 
bands in a unit cell 2a and the second term is related to the intraband transition, e is the electric charge and a 
are the Pauli matrixes. We neglect the relative distortion 77 (= 2u/a) in the dipole operator because it is relatively 
small in the optical contribution. 

Due to the fact that tt electrons in the SSH model are non-localized, the dipole approximation [^,^ is no longer 
valid in the extended states and will lead to wrong results as pointed out by some authors. Q The Fourier transform 
of Eq. (Q) to coordinate space shows that the transition dipole is related to the electron hopping to all the other sites 
besides the nearest neighbor sites. Thus the dipole approximation fails for the extended states in periodic systems. 
Because of the failure of dipole approximation through the polarization operator P, computations show a magnitude 
difference of 10^ in x^""^-* and 10^ in x'"^^ [0 and quite different shape in spectrum compared with the results through 
the dipole operator D and the experimental values in trans-(CH)x, although the position of some resonant 

peaks may be correctly obtained. 

LO response by D-D: The LO susceptibility x^ssni^^'^i) obtained from Eq. (||) and Eq. (|^): 



T r< i 



d_ 

"dk 



G{k,Lu)t-^[G{k,uj-uji)] 



+ [p{k)G{k, Lo)a2G{k, uj~uji)]+ l3{k)'^2G{k, u;)(3{k)a2G{k, co-loA^, (5) 



where the Green function G(k,uj) = — ^ ^ 



By Eq. (|), we have XssHi^)= Xslni-^^^)- 



,(1) 
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with ujk = e{k)/h and e = 0^ 



e\2taa) 
27rA2 



where x = huJk/A, z = 7ia;/(2A) and 5 = A/{2to). 

If the continuum limitation is applied, that is, S - 
XTLAii'^) under the TLM model [H as follows: 



dx 



[{1 - S^x^){x^ - l)]ix^{x^ - zJ)' 



and 2^00 — > hvp-, the above intergral gives the LO susceptibility 



where 



e^hvp 
"27rA202 



(6) 



arcsin(z) 



cosh^"'"(z) 



(z2 < 1), 
(z2 > 1). 



(7) 



The conductivity (tIuj) given by ~iuj ■ x^^^ , is exactly the same as based on J- J correlation. |l4] However, we should 



point out that the direct computation based on the static current from J -J correlation under TLM model |2yi4[| shows 
no first term in Eq. (^, [ pT| ZED in real part of J- J correlation is obvious although the correct imaginary part 
still can be given. These difficulties have never been clearly addressed previously, provided the reason that the static 
current Jq is still valid for obtaining the correct imaginary part in the LO absorption. To include the IFCs by changing 
k ^ k — eA/fi in the static current operator Jo(k), ZED could be solved to give the same result as that under D-D 
correlation. But the complicated way to include the lECs compared with simple D-D correlation already makes 
J-J correlation impractical even for the LO response. Fortunately, the lECs only contribute for the real part of J-J 
correlation and have no influence on the absorption. Attempts to obtain lECs directly from TLM Hamiltonian fail 
because it is forbidden to include term by the model. |^ In the discrete SSH Hamiltonian, we have the chance to 
include the IFCs through Peierls substitution, Q but a straightforward computation easily shows that they are not 
correct IFCs to cancel the ZED. It gives us an impression of the difficulty to obtain the correct IFCs besides the 
complexity of handling their contributions to the optical response even in well-defined 1-d periodic models. From the 
above examples, the feasibility of J-J correlation in more general models will be questioned and the application of 
D-D correlation is more reasonable under approximate models to obtain the correct results. 
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NLO susceptibilities of trans-(CH)x chain: There are many elegant works in discussing NLO susceptibilities of 
polymer chains. | |ll| , p^| -pOt The TP obtained from J- J correlation was doubted in the literature p5|-p0| since it 
is forbidden by momentum conservation and parity consideration in both TLM and SSH models. Based on the D-D 
correlation, our analytical results of THG show explicitly that TP no longer exists in both the SSH and the TLM 
models. This unphysical TP is caused by the same reason - the improper use of the static current Jo and omission 
of the IPCs in periodic chain. Although the reason for J-J correlation is very clear, we only give the computational 
results based on D-D correlation. It is expected that if the correct IPCs are considered, J-J correlation will give 
the same results as D-D correlation although much more complicated computations are inevitable. After a similar 
definition as Eq. (||) and tedious derivations, the new result of THG per unit length under SSH model for infinite 
chains is recovered as: 



(3) 

Xo 



45 
128 



dx 



47 - 48(1 + (52)a;2 + m'^x^ _^ 3(1 - 5'^x'^){x'^ 



1) 



+ 



8a;8(x2 - ^2) 

9 [47-48(l + ,52)a:2 +4 
8x8 (x2 _ {2,zY) 



x^(a;2 — z^Y' 

63(1 - (5^a:^)(a;^ - 1) 
a;6(a;2- (3z)2)2 



(8) 



where Xg 
oriented. 



(3) _ 8 e^no {ItQof 



hq is the number of chains in unit cross area, the polymer chains are assumed to be 

(1) 



45 TT A 

; and S are defined the same as in Xssiii^) 



Eq. (|^) is an elliptical integration and can be numerically integrated if one change x 



_ : + ie in considering the 

life-time of the state. ]lj,|l^ Por polyacetylene, by choosing to = 2.5eV, A = 0.9eV, uq = 3.2 x 10^"* cm^^, a = 1.221 
and e ~ 0.03, we have 6 = 0.18 and Xo^^ « 1-0 x 10"^° esu. The absolute value of xlsni^) is plotted in Fig. 1. 



It shows the good agreement with the experimental value 
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L5| around z = 1/3. 
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FIG. 1. \Xssh{^)\ for e=0.03 with z = hLu/{2A). 
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FIG. 2. Computed D-D value (solid line) vs. Jo- Jo value (dashed line)of Ixtlm (f^)! with z = huj/{2A). 
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Let (5^0"'" and e ^ 0+ in Eq. (p[), we obtain the analytical result of THG under TLM model as follows: 



^THG ^(3) 45 14 4 (37- 24z^) {I - 8z^) 1 



where /(z) and Xo^^ defined in Eq. (g) and Eq. (|). The comparison between our result(£'-I?) and Wu's result( Jq- Jq) 
[ pl| on absolute value of X^lf^ plotted in Fig. 2. 

The TP disappears in our analytical results, which is more reasonable for the physical situation and consistent 
with the previous numerial computations. |16-20| It also implies that the e-e interactions, p9| , p0| disorders, quantum 



fluctuations or finite chain size effects | |17[ should be taken into account to explain this experimentally observed TP. 

Fig.l shows another new resonant peak z = 1 at ratio of 1/10 of z = 1/3. This new peak hasn't been reported by 
the experiments ||l^ because of the experimental scanning range. The cancellation of ZFD by Jq- Jo correlation [|llj is 

(3) 

actually a coincidence under TLM model, because in the SSH model, we find ZFD in Xssh through Jq- Jo correlation. 
pl| Keldysh formalism |l^ is not necessary to apply in this equilibrium system. |l^ From THG of trans-(CH)2,, D-D 
correlation is superior to J- J correlation to obtained NLO susceptibilities. 

As a conclusion, the principle about the equivalence oi D-D vs. J- J correlation is still correct under the approximate 
models while practicality favors D-D over J- J correlation. Although our discussion is chiefly based on l-d polymer 
chains, the main results in this letter can be generalized for 2-d or 3-d periodic systems. 
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